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DISORDER IN CRYSTALS [m 


SYNOPSIS 


A general survey of order - disorder phenomena in 
crystals is presented. Several types of disorder are consi- 
dered, with special emphasis on disorder of position, 
disorder of orientation and magnetic disorder. A review 
of the main experimental techniques used in these stu- 
dies is given, as well as a general discussion of some 
theoretical treatments. 


A. TYPES OF DISORDER 


Any crystal, however perfect, is in a sense always 
disordered. It can only be regarded as achieving zero en- 
tropy at OK, the absolute zero, a temperature which 
(according to the third law of thermodynamics) can 
never in fact be reached. However, the non-zero entropy, 
and hence the disorder, which a crystal inevitably pos- 
sesses is in general a result of the distribution of its 
energy among the various possible vibrational and other 
energy levels available to the particles in the lattice. In 
this article we shall consider more specific kinds of 
disorder, namely the following: (1) Disorder of position. 
(2) Disorder of orientation. (3) Magnetic disorder. 


(1) Disorder of position. This can arise in more 
than one way. 


(a) Two or more kinds of particle can be ran- 
domly arranged on a given set of lattice sites. This situa- 
tion is encountered, for example, in alloys. 1.2 Thus, the 
alloys CuZn and FeAl exist in an ordered low tempera- 
ture form, and a high-temperature disordered form. In 
each of these two alloys, the lattice sites are the same in 
both ordered and disordered forms, namely those of a 
simple body-centred cubic lattice. In the ordered form of, 
say, FeAl, each Fe atom is surrounded by eight Al atoms, 
and vice versa, whereas in the disordered form the two 
kinds of atom are mixed up on the available sites, so that 
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RESUMO 


Apresenta-se uma revisão do problema da ordem- 
-desordem em cristais, focando-se em especial os casos 
de desordem de posição, desordem de orientação e 
desordem magnética. Faz-se uma discussão dos métodos 
experimentais mais utilizados no estudo destes proble- 
mas, bem como dos tratamentos teóricos mais correntes. 


an atom of either kind is, on average, surrounded by four 
Fe atoms and four Al atoms. In certain alloys, only some 
of the atoms are involved in the disorder. For example, 
in the alloy Fe, Al two thirds of the Fe atoms have fixed 
positions in the lattice in both the ordered and dis- 
ordered forms, while one-third of the Fe atoms and all of 
the Al atoms are mixed up on the remaining sites in the 
disordered high-temperature form. 


(b) There are more sites in the lattice than particles 
available to fill them. A famous case of this is ice, where 
the number of possible sites for the protons is twice the 
number of protons. Another example is the solid 
Hg(NH,), Cl,, which has a cubic unit cell3. Here, the 
Cl-ions have fixed positions (at cube corners), and the N 
atoms have fixed positions (at the cube centres). The Hg 
atoms are situated at the face-centres of the cubes, but 
there are six times as many of these positions as there are 
Hg atoms to fill them. The Hg atoms are randomtly distri- 
buted among these positions, so that the chance of any 
particular face-centre being occupied is 1/6. An extreme 
example of positional disorder is encountered in silver 
iodide, Agl.4 At about 146ºC this undergoes a transition 
into a high temperature form with a remarkably high 
electrical conductivity, (comparable with that of a fused 
salt.) This is because in the high-temperature form, 
whereas the iodide ions occupy fixed positions in the lat- 
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tice, the silver ions move freely between them. One might 
describe the high-temperature disordered crystal as fluid 
silver ions moving in a solid formed by the iodide ions. 


(c) Although in this article we shall not explicitly 
consider lattice defects and non-stoichiometric solids, it 
is worth noting that thermal agitation can create vacant 
lattice sites in a crystal (so-called Schottky holes or 
defects), and that these can make it possible for the 
particles in the lattice to move from one site to another. 
More information on this will be found in Section F. 


(2) Disorder of orientation. This arises when a 
diatomic or polyatomic molecule or ion has more than 
one distinguishable orientation in the crystal lattice. (The 
possible orientations must be recognizably different with 
reference, say, to the crystal axes.) Examples of this are 
very common and are to be found in simple molecular 
crystals (e.g. HCl, H,S, N,, O, CH,), in more com- 
plex organic crystals of which examples will be given 
later (Sections E and F), and in ionic solids (e. g. certain 
salts containing ions such as NH+ « ENT. NO” a 
A simple illustration is provided by ammonium choride. 
The structure of this, in both ordered and disordered 
forms, is like that of CsCl. In the tetrahedral NH* ion, 
the positive charge is distributed among the hydrogen 
atoms, and it is natural to expect that each N — H bond 
will point towards a Cl” ion. Any one NH : ion, situa- 
ted at the centre of a cube of eight Cl- ions, can so 
direct its four N — H bonds in just two ways (Fig. 1). 


Fig. 1. Showing the two orientations available to an NHY ion 


at the centre of a cube of Ci” ions in NH CL. 


In the low-temperature ordered form, all NH + ions 
adopt the same orientation, that is they are all parallel. 
(This is because the electrostatic interaction energy bet- 
ween adjacent NH o ions is lower in this arrangement). 
In the high-temperature disordered form, each ion makes 
use of the two possible orientations available to it, so 
that at any instant two adjacent tetrahedra may be 
parallel or they may be antiparallel. 

As another example of orientational disorder we 
may cite the case of deuterium bromide?, DBr which 
exists in three forms: | (the high-temperature form), Il 
(the middle form), and Ill (the low-temperature form). 
[This, incidentally, illustrates the more or less accepted 
way of numbering different phases]. In DBr-lll, all atoms, 
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both D and Br, occupy definite positions in parallel 
zig-zag chains, each chain being of the kind 


In DBr-ll, each DBr molecule has two possible orienta- 
tions, and can direct its D atom towards one of two 
neighbouring Br atoms. DBr-l is a face-centred cubic 
lattice in which any one molecule has twelve nearest 
neighbours, and it may direct its D atom towards any one 
of these, 

The limit of orientational disorder is free rotation 
of the molecule or ion. This was at one time thought to 
be common, and transitions at which free rotation was 
believed to begin in a lattice were often spoken of (and 
sometimes still are) as 'rotational transitions'. It now 
seems fairly certain that free rotation, in the literal sense, 
is very rare. One of the few cases in which it actually 
happens is in a crystal of n-H, , (that is, normal hydro- 
gen consisting of 3/4 ortho-molecules and 1/4 para- 
-molecules), where, near the melting-point at least, the 
ortho-molecules can undoubtedly rotate freely in the 
lattice. 


(3) Magnetic disorder. lons of transition metals 
and rare-earth elements which have unpaired electrons 
behave as elementary magnets, and if a lattice contains 
such ions these must interact to produce, at sufficiently 
low temperatures, an ordered spin system. The ordering 
can be of two basic kinds. Either the spins can line up in 
the same direction, in which case the magnetic moments 
reinforce each other to produce a ferromagnetically 
ordered system, or else the spins order themselves in 
a selí-cancelling way to give an antiferromagnetically 
ordered structure. A variation of the latter type of orde- 
ring is when the crystal contains two sorts of magnetic 
particle with different moments, as in Fe,0,, so that 
even if the spins of the two kinds of particle adopt an 
antiparallel alignment they do not completely cancel. A 
domain of cells therefore has a resultant magnetic mo- 
ment. Such a solid is termed a ferrimagnetic. The magne- 
tically disordered state in a transition or rare-earth com- 
pound is simply, of course, the normal paramagnetic 
condition, in which the spins are randomly orientated but 
tend to align themselves in an external magnetic field. 
The temperature at which they adopt one of the ordered 
arrangements just described can vary greatly. It may be 
several hundred degrees Kelvin, or it may be below 1 K. 
There is one very common molecule where the possibility 
of magnetic ordering in the crystal arises, namely oxygen. 


(4) The above three types of disorder do not ex- 
haust the possibilities. Sometimes the particles in the 
lattice (certain rare-earth ions are examples) have 
available to them two or more electronic energy levels, 
and as the temperature rises the ions gain access to 
the higher levels with absorption of energy. And some 
times interesting thermal phenomena are encountered at 
very low temperatures which have their origin in nuclear 
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spin and tunneling through energy barriers. Examples of 
this have been found in methanef and in nickel hexam- 
mine iodide.? 


(5) We may end this section with some general 
observations. One is that the distinction we have made 
between positional and orientational disorder is some- 
times rather arbitrary. For example, the disorder in ice 
can be thought of as a matter of positional disorder 
of the hydrogen atoms among the more numerous po- 
sitions which they can occupy in the lattice (subject to 
the condition that the unit in the lattice is H.O and 
that only one hydrogen atom can lie on the line between 
any two adjacent oxygen atoms), or alternatively as a 
question of the orientational disorder of the bent H, O 
molecule in the lattice. Anothed point worth making is 
that a crystal can sometimes simultaneously display more 
than one kind of disorder. Thus, in lithium iodide mono- 
hydrate (Lil, H, O), which has a perovskite structure (Fig. 
2), only the iodide ions and the centres of the water mole- 
cules have definite positions.8 The Lit ions are statisti- 
cally distributed among the face-centres of the unit cube 


(V) Q 
H,0 8 


Fig. 2. The structure of LilH O. The Li+ ions are situated at 

the face-centres. The chance that any one face-centre is occupied 

by an Li+ i10n is 1/3. The H O molecules are orientationally 
disordered. 


(and so are positionally disordered), while each water 
molecule (which is at the centre of the cube, surrounded 
by eight |” ions) is orientationally disordered, as there 
are various ways in which it can direct its two hydrogen 
atoms towards two of these eight |” ions. 

Finally, a distinction should be made between 
long - range order and short - range order in a crystal. In 
saying that at high temperatures the alloy FeAl is com- 
pletely disordered, we mean that on average any site 
in the lattice (picked at random) is equally likely to be 
occupied by an Fe atom as by an AÍ atom, and that 
there is no correlation between the nature of the partic- 
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ular atom at any one site and of those at sites some 
distance from it. (In the perfectly ordered crystal, by 
contrast, it is sufficient to say this particular site is 
occupied by an Fe atom to determine which kind of atom 
occupies every other site in the crystal, however remote). 
However, the nature of the low-temperature ordered 
crystal shows that an Fe atom prefers (on energy 
grounds) to have Al atoms rather than Fe atoms as its 
nearest neighbours. Hence, by the well-known Boltzmann 
principle, which requires that in any system of particles 
there must always be some preference for configurations 
or arrangements of lower rather than higher potential 
energy, even in the disordered alloy there will be some 
tendency, however slight, for an Fe atom to prefer Al 
atoms as its nearest neighbours. This gives rise to some 
local, or short-range order. In general in this article we 
are considering average, or long-range, order. 


B. TEMPERATURE AND DISORDER 


(1) If a crystal, disordered in one of the ways des: 
cribed above, is cooled down, it generally (though not 
always) becomes ordered, essentially because the order- 
ed condition is a state of lower energy. It is perhaps 
better to think of the process in reverse, that is as an or- 
dered crystal being heated up from a low temperature. 
The change to a disordered crystal can then occur in 
two ways. 

(a) Isothermally, in a so-called first-order phase 
change. An example is provided by the || -—» I phase 
change in HCl at 98.4 K. Here the structure changes from 
a face-centred orthorhombic lattice in which the mole- 
cules are orientationally ordered, to a face-centred cubic 
structure in which they have twelve-fold orientational 
disorder. At such a transition there is an abrupt change 
in the enthalpy H, entropy S, and volume V, and in the 
first differential coefficient, (cG/dT)p. of the Gibbs 
energy — hence the term 'first-order”. 

(b) Gradually, that is over a range of temperature. 
Such a transition is revealed by a so-called 'anomaly' in 
the heat capacity Cp . Over a range of temperature (the 
extent of which can vary greatly) Cp is higher than 
it would “normally" be. The “extra” heat capacity, plotted 
against temperature, looks something like the Greek 
letter 2, and transitions of this kind are frequently des- 
cribed as ?-type transitions. The order-disorder transition 
in ammonium choride which we have already mentioned 
is of this kind. Two >-type transitions are shown by pota- 
ssium cyanide?, KCN (Fig. 3). Such gradual transitions 
are very common, occurring in solids of almost all 
chemical types (and the one known transition in a purely 
liquid system, that between the two forms of liquid 
helium, is of this kind). >-typo transitions may, or may 
not, involve a change in the basic crystal structure. Thus, 
in ammonium chloride, this is simple cubic both above 
and below the transition. In KCN the two lower forms are 
orthorhombic, but the high-temperature form is face-cen- 
tred cubic. 

The essence of a gradual transition is that it is co- 
-operative in nature. This requires that the change 


175 


30 


> 
ja) 


C,. cal./mol. dee. 
A 


0 so 100 


order -» disorder becomes energetically easier, the more 
disorder there is in the system. The change may be 
thought of as analogous to a chemical reaction in which 
the “reactants' are the ordered crystal and the 'products' 
the disordered crystal, and in which the change from 
'Teactants' to 'products' is endothermic — but with this 
essential difference from a conventional chemical reac- 
tion, that in the change from order to disorder the heat of 
the 'reaction' decreases as the “reaction” proceeds, be- 
coming zero when it has gone over completely to the 
products" side. 

Order-disorder transitions in crystals cannot really 
be divided sharply into two kinds, namely either first- 
-order (wholly isothermal), or completely gradual. Pro- 
bably many transitions begin gradually but reach comple- 
tion isothermally. Whether this is so in any particular 
case is best examined experimentally by using a property 
which can be accurately measured at a series of fixed 
temperatures, for example volume or dielectric constant, 
rather than the heat capacity (the measurement of which 
means heating the sample over a range of temperature, 
even though this may be small). Thus, a careful dila- 
tometric study of the transition in NH ,ClI (often quoted 
as a typical »type transition) showed that the transition 
begins gradually but is completed isothermally.10 (By 
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Fig. 3. Plot of the molar heat capacity (Co ) of KCN against temperature. 


contrast, the corresponding transition in ND Cl is com- 
pletely gracdual!! — an example of the sensitivity of the 
details of transitions of this kind to the nature of the 
particles in the crystal.) 


(2) Sometimes a disordered crystal remains disor- 
dered when cooled down to the lowest possible tempera- 
ture. Examples of familiar substances which behave in 
ihis way are H, O, CO, NO, N,O and Na ,SO,, 10H,0. 
Such a crystal then appears to disobey the third taw 
of thermodynamics, in that it retains entropy at OK. We 
shall call this entropy Sres' (res = residual). Sres can 
be determined by comparing the value of Scal” the appa- 
rent calorimetric entropy of the substance in some defi- 
ned state (e.g. at 25C and 1 atmosphere pressure) 
with the 'true' value of the entropy in this state. This 
latter quantity will be given for a molecular substance 
like CO by Sspec the entropy calculated from spectros- 
copic data using the standard methods of statistical me- 
chanics, or alternatively for a substance which is solid or 
liquid in the defined state by Seq the value obtained 
from some suitable equilibrium study. Sçaj Will be given 


PT 
by | Cp Ir dT Ee >4H/T where Cp is the experimentally 
J O 


determined heat capacity and 24H /7 takes account of 
the entropy increases at phase changes or transitions. 
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Examples of the determination of Seq will be given later. 
Sres is, of course, the difference between S spec (or 
A ) and Scal 

f a crystal remains disordered down to OK, it is 
generally not because the ordered and disordered con- 
figurations have equal energies. The ordered state will 
still be expected to have the lower energy, but the energy 
difference is now relatively small, so that it only becomes 
significant at relatively low temperatures — at temperatu- 
res so low, in fact, that by the time they are reached the 
particles in the lattice find it impossible to undergo 
the necessary rearrangement to give the ordered crystal 
in a reasonable period of time. So a crystal like solid 
NO, which is still disordered at OK, is not really in a 
true equilibrium state. The disorder has been frozen in, 
as it were. 


C. EXPERIMENTAL METHODS 


A wide variety of experimental methods is now used 
in investigating disorder in crystals. The most important 
of these are briefly listed below, with examples of the 
sort of information they provide. 


(1) Thermodynamic studies. These are inevitable at 
some stage of the investigation, and indeed most order- 
-disorder transitions have been discovered by heat ca- 
pacity investigations, which have either revealed a 
first-order phase change or a gradual >-type transition. 
If the transition is gradual, the heat capacity cetermina- 
tions can lead to information on how the extra entropy is 
acquired. Thus, for magnetic transitions in particular, 
sometimes only a relatively small fraction of the total 
entropy taken in during the transition has been gained 
at the temperature at which the maximum in the heat 
capacity is reached, implying that considerable order still 
persists above this temperature (Section H). Moreover, 
as we shall see, entropies of order-disorder transitions 
are sometimes (though by no means always) capable of 
a simple and informative statistical interpretation, And, 
as already mentioned, the study of crystals which remain 
disordered down to OK requires the experimental ther- 
modynamic work necessary to give Scal' and (if appro- 
priate) Sea 

In experimental thermodynamics, pressure is much 
less often varied than temperature, and undoubtedly a 
good deal of interesting work remains to be done on the 
effects of high pressure on disordered crystals. One sys- 
tem which has been very thoroughly investigated under 
pressure is ice, and we shall later refer briefly to some 
of the high-pressure forms of this solid. 


(2) Diffraction methods. Since we are dealing with 
solids, X-ray diffraction methods are of great importance. 
They have their limitations, however. They do not, for 
example, in general distinguish between the genuine free 
rotation of an ion or molecule (though on other evidence 
this is almost certainly a very rare occurrence, as we have 
said) and its orientational disorder between a limited, 
definite number of directions. In seeking a thorough un- 
derstanding of particular cases of disorder, it is obviously 
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sensible and desirable to study crystals in which the ions 
or molecules are relatively small. Solids composed of 
small molecules are often gaseous at ordinary tempera- 
tures (e.g. N,,0O,,F,), so that the experimental di- 
ffraction work on such crystals must be carried out at 
a low temperature. An example of what can be achieved 
is shown in Fig. 4, which represents the crystal structure 


Fig. 4. The structure of the high-temperature form of fluorine, 
tr, = = 8=-F lh The balls represent F. molecules which 
have orientational disorder in three dimensions. The remaining 
F, molecules are orientationally disordered in two dimensions, 
- (i.e. in a plane). 


of the high-temperature disordered form of fluorine, as 
determined by Lipscomb and co-workers!2, This is a 
particularly interesting structure, as the molecules fall 
into two groups, those of one group being orientationally 
disordered in three dimensions (spherically disordered), 
those of the other group in two dimensions only (cylin- 
drically disordered). The limitations of X-ray diffraction 
in locating light atoms in a solid, particularly hydrogen 
atoms, are of course well-known, and in this respect 
electron diffraction can be a somewhat superior techni- 
que. Thus, the disordered structure of ice has been con- 
firmed by an electron diffraction study.!3 But still more 
valuable for locating hydrogen atoms is the technique of 
neutron diffraction, usually performed on the deuterated 
compound. Elegant work has been carried out, for exam- 
ple, on hydrogen and deuterium halides by Sandor and 
his colleagues, in some case with a degree of refinement 
which actually makes it possible to eliminate the possi- 
bility of free rotation in favour of reorientation between 
a definite number of orientations. Thus, Sandor and 
Farrow!4 have shown that in the high-temperature face- 
-centred cubic form of DCI, the molecules are not freely 


177 


rotating but are distributed among twelve possible orien- 
tations (each molecule having twelve nearest neigh- 
bours), (Fig. 5). Furthermore, since the neutron has a 
magnetic moment and interacts with an unpaired elec- 
tron, neutron diffraction is a powerful tool in investiga- 
ting magnetic order in crystals, which has been much 
used in studying compounds of transition metal and 
rare-earth metais. With solid oxygen, which exists in 
three forms, there is the possibility of both magnetic 
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Fig. 5. Neutron diffraction study of DCI — |, the cubic high-tem- 


nerature form. |F| is the intensity of the diffraction maxima, and 4 

is the diffraction angle. A shows the experimental results, B is 

the calculated pattern for a disordered model in which each 

DCi molecule has twelve possible orientations, and C is the 

calculated pattern for randomly orientated (freely rotating) mo- 

lecules. It will be seen that the agreement with A is better for B 
than for C. 
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and orientational disorder. A neutron diffraction study of 
the lowest temperature form Il ( =z ) and of the 
middle form ll ( = 5 ) made it clear that Ill, but not II, 
is magnetically ordered,1º (Fig. 6). It will be seen from 
Fig. 6 that the two low-angle peaks shown by phase Ill are 
not given by Il, and their appearance for Ill must be due 
to the magnetic ordering of the oxygen molecules in this 
phase. It will be noticed, however, that for phase Il the 
experimental curve for the intensity of the diffracted 
neutron beam shows a slight bump in place of these two 
peaks, implying that there is still a tendency towards 
local magnetic order in this phase, though, of course, no 
long-range order. 


(3) NMR (Nuclear Magnetic Resonance). This tech- 
nique has been especially valuable in giving information 
about the kinetics of the molecular movements associated 
with disordering processes, that is about the rate at 
which molecules or ions can change their orientations or 
move from one lattice site to another, and about the 
activation energies for these processes.16 Use has been 
made of the effect of temperature on line-widths, second 
moments, and spin-lattice (T, ) and spin-spin (T, ) re- 
laxation times. As just one example of the sort of 
behaviour encountered on passing from an ordered to a 
disordered crystal, we reproduce in Fig.7 the temperature 
dependence of the linewidth of  pentaerythritol, 
C(CH, 0H), .17 This substance undergoes a first-order 
transition to a disordered phase which exemplifies a type 
of solid called a plastic crystal (Section F). The line- 
width starts to decrease below the transition-temperature 
due to the gain by the molecules of orientational freedom. 
The activation energy for the reorientation process is 25.2 
kcal mol-1. In the plastic phase the linewidth decreases 
still further to relatively very small values. This further 
decrease is ascribed to the ability of the molecules to 
move rapidly from one lattice site to another, with an 
activation energy estimated to be 24.3 kcal mol-?, 


(4) Dielectric constant studies. For molecular crys- 
tals with polar molecules (e.g. HCl, H, S, H, Se), the 
acquisition of orientational freedom can often be revealed 
by studying the temperature-dependence of the dielectric 
constant (d.e.c.)18, For the ordinary substance with 
polar molecules which exists in just one form (e. g. ni- 
tromethane, CH, NO.), the d.e.c. of the crystal is 
relatively small, since the molecules in the lattice are 
unable to change their orientation when placed in an 
external electrical field. The d. e. c. of such a substance 
rises abruptly at the melting-point. For HCl, on the other 
hand, the d.e.c. increases sharply at the first-order 
transition from the orientationally ordered form Il to 
the orientationally disordered form |. When | melts, the 
d. e. c. actually drops cue to the expansion on melting. 

Valuable information on the kinetics of the reorien- 
tation of polar molecules in crystals can be obtained 
from studies of the phenomenon of dielectric dispersion 
(essentially the change of d. e. c. with frequency). Exam- 
ples of this will be found in a paper by Cole and Ha- 
vriliak.19 


(5) Infra-red Raman spectra. The infra-red spe- 
ctrum of a crystal can be studied both in the near infra- 


TÉCNICA 414 


Reinos rue p 


u “ds 
o. 
a |) 
$ 
* 


Aspectos dos 2 guindastes pura trabalhos com gancho e 
colher com a capacidade de carga de 6 toneladas e 21 me- 
tros de alcance em serviço no caís da SAPEC em Setúbal. 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 


TURBINAS HIDRÁULICAS ——— Fabrico segundo licença de A.C.M. de Vevey, S.4. 
TURBINAS A VAPOR ————.— —— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR ——————— Projecto e fabrico segundo licença de Foster 


Wheeler, Co. 
EQUIPAMENTOS E INSTALAÇÕES e 


INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGUE :..-. 


ALVERÇCA DO RIBATEJO- PORTUGAL 
| TECNICA XII 


. BY 
EM: 
Ro 
| Wa : 


Maias a 
A utilização de betão pronto impóôs-se 
ao dono da obra, 
ao técnico responsável 
e ao empreiteiro 
porque garante qualidade, 
rapidez e economia 
em relação 
a um betão fabricado em estaleiro. 


Consulte-nos. 


Os nossos serviços técnicos 

e comerciais 

terão o maior prazer em prestar 
todo o apoio 

e colaboração 

que nos forem solicitados. 


UMA NOVA MARCA 
DE BETAO PRONTO 


JOIVIA NET: 


TECNICA XIV 


ES SP SUIÇO-PORTUGUÊS | 286/439/AR 


VÁLVULAS 


de cunha, 
retenção, regulação 
e passagem 


Segundo normas DIN e ASA, 
para líquidos, vapores e gases 
(incluindo Cl, NHs e SO» ). 
Execução em ferro fundido, aço 
vazado, aço inoxidável 

e ligas especiais. 

Diâmetros nominais de 10 a 600mm. 
Pressões nominais de 4 a 40 kg./cma. 


ESTUDOS +» PROJECTOS 
ORCAMENTOS 
ASSISTÊNCIA TÉCNICA 


SUCIEDADE PORTUGUESA 
DE VÁLVULAS, 1º 


Rua da Academia das Ciências, n.º 5, 
Lisboa-2 = Telefone P.P.C.A. 320011. 


LS 


15000!- "a Eai 


Magnetic lines 


= 
8 


Scattered neutron counts 


5000| - 


Scatterinq angle 26 


Fig. 6. Neutron powder patterns from q — CG (=Ill)at 4.2K and from 8-0 (=|) at 27K. The first 

two lines from q — O | are not observed on X — ray powder photographs, and indicate the presence of 
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-red (the region where intra-molecular vibration frequen- 
cies characteristically reveal themselves) and in the far 
infra-red. The latter region is being increasingly used 
now that the technical diffilculties have been largely over- 
come, and is particularly informative since it provides a 
means of getting direct information on the frequencies 
with which molecules rock (that is, undergo torsional 
oscillations) in the lattice, Generally speaking, both the 
near and far infra-red spectra of an ordered crystal show 
more detail or structure than that of the same substance 
in a disordered state. This is because the environment of 
a molecule or ion in the ordered crystal will be more re- 
gular or uniform than that for a disordered crystal, in 
which the environment will tend to vary at any instant 
from one molecule to another, so that the observed spec- 
trum will represent an average state of affairs. This point 
is illustrated by Fig. 8, which are shown the far infra- 
-red spectra of the three forms of DBr20, already men- 
tioned in Section A(2). The spectrum of the ordered 
form Ill reveals both lattice vibration frequencies (the 
two peaks at the lowest frequencies) and torsional 
oscillational frequencies (the two central peaks). In the 
partially disordered form Il, there is less structure, but 
the lattice vibration frequencies are still visible, while 
the still more disordered form | has a much more diffuse 
spectrum, rather like that of a liquid. In Fig. 9, part of 
the near infra-red spectrum of one of the high-pressure 
forms of ice, namely Il, is compared with that of form 
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| (the form stable at ordinary pressures) and forms Ill 
and V.21 In |, Ill and V the protons are positionally 
disordered, whereas in |l they are ordered. Again, it will 
be seen that the spectrum whith the most structure is 
that of the ordered form ll. 


Rather less use has been made of Raman spectra, 
but the spectrum of crystalline hydrogen within a few 
degrees of its melting-point supplied an elegant proof of 
the statement made earlier (Section A (2)) that this is a 
rare example of a lattice in which molecules can really 
rotate freely.22 The Raman spectrum of para-hydrogen 
(for the molecules of which the rotational quantum 
number, J, = 0,2,4, ...) shows a line with A = 355.5 


3 
cm-1, while that of normal hydrogen (consisting of p 


ortho-hydrogen and E para-hydrogen) gives this line and 
one with A = 589 cm:l. These two frequencies are 
almost exactly those for the rotational transitions 
J=0--> 2andJ=1-+3 respectively, being virtually 
identical with the corresponding Raman frequencies for 
the molecules in the gaseous state. 


(6) Magnetic susceptibility studies. In addition to 
the neutron diffraction studies already mentioned, the 
temperature dependence of the magnetic susceptibility x 
can give valuable information on the type of magnetic 
ordering ultimately adopted on cooling a paramagnetic 
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Fig. 8. The far infra-red spectra of solid DBr. A and B, phase Ill at 85K. C and D, phase Il at 95K. 
E, phase | at 125K. 


180 


TÉCNICA 414 


á A, | A 
DARE 


ep 
2 
S 
» 4 
sr € 9 é 
2 
4 
d h 
3 a 
e 
] | 1 , 1 , ' - 1 
3400 3300 1 2500 2400 
v km 


Fig. 9. “oH (HDO) and * (HDO) in various phases of ice as 


ob 
mulis in isopentane at — 100ºK. The ices used were: Curve a; 
Ice lh, 95% D O 5% H O: Curve b; Ice |1,95% D O 5% H O: 
2 q 2 
Curve c; lce III, 95% DO 5% H,0: Curve d; Ice V, 95% D,O 
5% H, O. Curve e; Ice lh, 95% H O 5% D, O. Curve f upper; 
Ico Il, 95% H,O 5% D 90: Curvo i lower; Ice Il, 99% H, O 1% 
D. O: Curve g upper; Ice Ill, 95% H, O 5% D.O: Curve g lower; 
Ice II, 99% H ,O 1% D,O: Curve h; Ice V, 95% H, O 5% D,O. 
The lower curve in f is plotted 0.05 units lower than observed, 
and curve h is plotted 0.1 units higher than observed, 


crystal. The characteristic dependence of x on T, de- 
pending on whether the solid becomes ferromagnetically 
or antiferromagnetically ordered, is shown in Fig. 10. 


(7) Slow neutron scattering studies. This is an 
interesting technique, likely to be put to increasing use 
in studying disordered crystals. Neutrons from an atomic 
pile are slowed down until their energies are comparable 
with the energies associated with the vibrations and 
torsional oscillations of the molecules in the lattice. 
Changes in the energy of the neutrons in the beam on 
passage through the crystal can therefore give infor- 
mation on the movement of the molecules. An example 
of this will be found in recent work on methane. 23 


(8) Radiotracer techniques. Some use has been 
made of these for obtaining information on the kinetics 
of the movement of molecules from site to site in a 
lattice, as for example in the high-temperature disorde- 
red phase of cyclohexane (the so-called plastic crysta- 
line phase).24 
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Fig. 10. Illustrating the characteristic dependence on temperature 

T of the magnetic susceptibility x of a paramagnetic solid (1), 

a solid which becomes ferromagnetically ordered (2), and a solid 
which becomes antiferromagnetically ordered (3). 


D. THEORETICAL ASPECTS 


In a general review of this kind it is impossible 
to deal adequately with the very considerable amount of 
theoretical work carried out on order-disorder problems 
in crystals, but attention will briefly be drawn to some 
of the aspects of the subject which theoretical treatments 
must consider. The equilibrium condition of a crystal is 
determined by the balance between two opposing fac- 
tors. One of these is thermal agitation, whith tends to 
produce a disordered state, The other is the interaction 
between the particles in the lattice, which tends to order 
the particles in a way which will minimize the energy 
of the system. The theoretical investigation of the ba- 
lance reached at a given temperature between these 
opposing factors is a matter for statistical mechanics. 
In general, different arrangements or different mutual 
orientations of the particles in the lattice will mean dif- 
ferent energies, and to deal with these some assumption 
has to be made about the way the particles interact, 
that is to say some “model” has to be adopted. Almost 
always the model is a simplification of the true state of 
affairs. One of the best known is the Ising model, origin- 
ally introduced to deal with magnetic interactions, but 
now used in a more general way (ses, for examples refe- 
rences 25 and 26. This at first sight seems a very simple 
kind of interaction, but even so its use can lead to mathe- 
matics of great complexity.) Sometimes the energy diffe- 
rences are small enough for them to be disregarded to a 
first approximation, but this can still leave 'configura- 
tional problems”, for example problems of the number of 
possible distinguishable ways of distributing particles 
on a lattice. These, too, can be of considerable difficulty. 
A famous calculation by Pauling27 of the number of 
acceptable ways W of distributing the 2N protons in 
a mole of ordinary ice gave the answer W = (3/2)N, 
But it is now known that this is not quite exact, though 
the value of kIn(3/2)N, = Rin(3/2) which it gives for 
the residual molar entropy of ice is believed to be correct 
to within about one per cent. (An exact value for W has 
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not yet been obtained, but a very good estimate is con- 
sidered to have been reached by computation)28. An- 
other configurational problem is the so-called dimer pro- 
blem. Imagine a primitive cubic lattice, which is to be 
occupied in a disordered way with dimers, that is rod- 
-like particles the two ends of which occupy two adja- 
cent cube corners. In how many ways can the dimers be 
fitted into the lattice? This problem has also so far resis- 
ted exact solution. 

With regard to the forces responsible for the orde- 
ring of a lattice at sufficiently low temperatures, the na- 
ture of these must naturally change from one type of 
crystal to another. Thus, that the FeAl alloy becomes 
ordered in the way described in Section A(1) must mean 
that Fe-Al interactions are energetically preferred com- 
pared with Fe-Fe and AI-Al interactions. It might be 
supposed that orientational order in a simple polar 
diatomic substance like HCl or DBr is due to dipole- 
-dipole interaction, and this may be so, but it is interes- 
ting that in the low-temperature ordered forms of these 
halides the adjacent molecules in the chain (see Section 
A(2)) are at about 90º to each other, which is not the 
arrangement which gives the strongest electrostatic 
attraction between dipoles. It implies rather that the 
attraction arises between the positive charge on the 
hydrogen atom and an unshared pair of p-electrons on 
the halogen atom, so that it may perhaps be regarded 
as a case of hydrogen bonding. In the orientationally 
ordered forms of some well-known diatomic substances 
(H,,N,.CO, F, ) there is little doubt that the ordering 
force is that between the molecular quadrupoles. The 
charge distribution is such that a molecule may be repres- 
ented thus: 


and an energetically preferred orientation between two 
such quadrupoles is the following: 


In their low-temperature ordered forms, all the four 
diatomic substances just mentioned have a cubic lattice 
made up of four sub-lattices, arranged so that this 
energetically favourable orientation can contribute to 
the stability of the lattice. It is significant that this struc- 
ture is the same as that of solid CO, , the molecule of 
which has a relatively large quadrupole moment. In solid 
tetrahydrides such as CH, , CD,, and SiH, , the syme- 
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try of the molecules is so high that the 'first' electrical 
moment the molecule possesses is an octupole moment, 
which is considered to be responsible for the orien- 
tationally ordered low-temperature forms of these subs- 
tances. 

In magnetically ordered solids, the forces producing 
the order are of yet another kind. The magnitude of 
these forces in, for example, transition metal oxides, 
sulphides and halides generally proves to be greater than 
ore might expect if the transition metal ions merely 
affected each other as if they were magnets interac- 
ting through space, and it is generally believed that in- 
teraction is transmitted through the electronic system 
of the intervenirg anions by a form of electron deloca- 
lization, or covalency, called 'superexchange'29, So one 
consequence of the study of magretic order-disorder 
transitions is that something may be learnt about the 
electronic structures of the solids concerned. 


E. ORIENTATIONAL DISORDER AND ENTROPY 


It is sometimes possible (though by no means 
always) to interpret the entropy gain associated with an 
order-disorder transition on the basis of the well-known 
equation: 

S = k Inw, 


where k is the Boltzmann constant and w is the number 
of microstates (sometimes called the number of com 
plexions) for the system.30, 31 Thus, if a diatomic or 
polyatomic ion or molecule goes from a condition in 
which it has n, distinguishable orientations to one 
in which it has a larger number n, of distinguishable 
orientations, the entropy gain per mole of the molecule 
or ion on this account (assuming that at any instance the 
ion or molecule can make full use of the orientations 
available to it) is k In(n,/n,)N, = R Inftn /n)). (The phra- 
se “on this account” is inserted because simultaneously 
with the gain in orientational disorder there may be a 
crystal structure change, or at least a volume change, as 
a result of which the vibrational frequency spectrum of 
the lattice alters, so that there is also a change in the 
vibrational entropy of the crystal. Generally speaking, 
such secondary contributions to the entropy change will 
be less, the less the difference between the crystal 
structures of the ordered and disordered forms). In the 
case of the » transition in ammonium chloride, in the 
low-temperature ordered form all the tetrahedral ammo- 
nium ions have one and the same orientation — they 
are all parallel to one another — so that n,= 1, 
whereas in the disordered form n, = 2, and hence the 
entropy gained on this account is R in2, = 1.38 cal 
mol-! K-1. This agrees with the experimental figure 
within the precision with which it can be measured.32 
The lower transition Il -—» Il in KCN (Fig. 3) involves 
no change in the orthorhombic crystal structure, and 
since the measured entropy gain of 1.32 cal mol! K-1 
is close to R In2, it is probable that each CN” ion in 
the low-temperature form Ill has a definite orientation 
(n,= 1), whereas in Il it has two possible orienta- 
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tions, CN— or NC”, (n, = 2), — sometimes referred 
to as 'head-to-tail' disorder. In the face-centred cubic 
lattice | where each CN” ion has six K* neighbours, 
it is likely that it has six possible orientations, so that 
the entropy gain at the Il — | transition due to the 
increase in the number of possible orientations for each 
anion from two to six is R In3, = 2.18 cal moll K1. 
The experimental value? is 2.01 + 0.03 cal mol! K-1, 

The substance benzthiophene33 has a >-type tran- 
sition just below room temperature, at which the entropy 
gain is 2.75 cal mol-1 K-1. This is almost exactly R In4, 
which implies that while below the transition all 
molecules have a unique orientation in the lattice, above 
the transition each molecule makes use of four possible 
orientations, as follows: 


Ss s 
Nevertheless, it must not be 
entropy changes of all order-disorder transitions, 
especially when they are first-order, can be so 
simply explained, Thus, the Il -—-» | phase change 


in HCl and DCI almost certainly takes the molecules 
from a lattice in which they each have a unique orien- 


12 


supposed that the 


10 


Number of transitions 
(Sp) co 


E 


Rin2 


tation to one in which they have twelve possible orien- 
tations. But the entropy gain at this transition is not 
R In12, = 4.94 cal mol! K-1, but only 2.89 cal mol! K-1, 
This 'discrepancy' is no doubt partly due to the fact 
that there is a crystal structure change and therefore 
some change in the contribution which the lattice vi- 
brations make to the total entropy. But in addition a 
molecule in the high-temperature form cannot choose 
which of its twelve possible orientations it will adopt 
independently of the orientations its neighbours may 
happen to have at that moment, since this would some- 
times lead to the arrangement Cl H... H— Cl, which 
is virtually impossible because of the repulsion between 
the hydrogen atoms. 

lt is characteristic of the structure of a salt that 
the cations and anions are arranged in alternate shells, so 
that if for example the anions are monatomic while the 
cations are diatomic or polyatomic, the latter are, as 
it were, protected to some extent from each other by 
an intervening shell of anions. This might be expected 
to reduce somewhat the sort of mutual interference 
which we have just illustrated by the example of HCHI, 
and therefore one might predict that the entropy gains 
associated with the transitions in salts involving orien- 
tational disorder would more nearly be given by R Inx, 
where x is a relatively small integer or a simple fraction. 
This is indeed found.34 Fig. 11 is essentially a represen- 


Rin3 Rin4 


AS 


Fig. 11. Distribution of the values of the entropy of transition AS for thirty-five transitions in salts for 
which AS lies between 0.8 and 3 cal mol"! K-!, The width of each colum is 0.2 cal mol! KI 
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tation of the distribution of 35 experimental values of 
transition entropies in salts. The grouping of these values 
round R In2 and R In3, and to lesser extent round R In4, 
is striking. 

When a crystal remains in a disordered state down 
to OK, the value of the residual entropy Sres can often 
be interpreted on the basis of the expression k Inw. 
Crystalline N, O is an example. 35 The molecules are 
linear, and apparently disordered as between the two 
possible orientations NNO and ONN. If the disorder 
were complete we should expect that Sçes would be 
R In2, = 1.38 cal mol! K-1. In fact, the experimental 
value (as given Dy Ss ecT Scal) is a little less than 
this, being 1.14 cal mol-! K-1, implying that there is 
some degree of local order, that is, that there is a slight 
preference for the neighbours of a given molecule to 
have one of the two possible orientations rather than 
the other. Crystalline nitric oxide, NO, presents an inte- 
resting case of residual disorder.36 The crystal is dia- 
magnetic, unlike the liquid and gas, and the unit in the 
lattice is a dimer, N.,O, . The mole of NO provides 
N/2 dimers, and the assumption that each dimer has 
two possible orientations leads to the prediction that Sres 
should be k In2N/2, = R/2 In 2, = 0.69 cal mol1 K“1. 


The experimental value, i.e. S — S is 0.72 


cal mol! K-1. 


spec cal” 


F. PLASTIC CRYSTALS 


This name has been given to a type of disordered 
molecular crystal (usually, though not always, of an 
organic compound) with relatively small entropies of 
fusion (usually less than 5 cal mol! K-1).37 Table 1 
gives a few examples of substances which form plastic 
crystals. (Some of the simple substances we have 
already mentioned, for example N, ,F, and CH, , could 
also legitimately be described in their disordered phases 
as plastic crystals). A more complete list has been 
given by Timmermanns38, who first recognized the exis- 
tence of this interesting type of crystal. He noted that 
the molecules forming them tended to be approximately 
spherical, including for example, tetrahedral molecules 
like neopentane (C(CH,), ), cyclic molecules like cy- 
clohexane, or cage molecules like camphor, and he pro- 
posed the term “globular molecules' to describe them. 
The plastic crystal is formed from a low-temperature 
phase at a transition which may be first-order or 2.-type 


TABLE 1 


Some examples of plastic crystals, showing: (1) the temperature at which the plastic phase is formed from the 
'“ordinary” low-temperature phase: (2) the melting-point of the plastic phase; (3) the corresponding entropies of traa- 
sition and fusion, (all 2. S values in cal mol-1 K-1); (4) for the first four examples, the nature of the transition. 
(F = first-order transition, 2» = gradual transition, for which the transition temperature is that at which the heat 


capaciy is a maximum). 


Carbon tetra- Ccil, 
chloride 
Carbon C(CH, ), 
tetramethyl 
Pentaerythrito!l C(CH, OH), 
Cyclohexanol C, H,, OH 
Camphor CH, 
CH E co 
HeCeH, 
CH;—C-CH, 
h 
Tungsten WF, 


hexafluoride 
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DE gia — DT :| «e 


"Ordinary" phase 


- LIQUID 


Plastic” phase 


225.5K(F) 250.3K, 

LS = 4.86 LS = 2.40 
140K( 2.) 256.5K, 

AS = 4,39 LS = 3.03 
457K(F) 539K, 

LS = 22.8 AS = 3.16 
265.5K(F) 299K, 

AS = 7.95 LS = 1.4 
250K, 453K, 

AS = 7.6 AS = 2.8 
265K, 275K, 

AS = 7.8 AS = 3.56 
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and at which there is a considerable entropy increase, 
greater than that which accompanies the melting of the 
plastic crystal. The difference between A S(transition) 
and / S(fusion) for pentaerythritol, (Table 1) is a 
striking example of this. Sometimes the entropy of 
fusion of the plastic crystal is surprisingly small; that 
of cyclopentane is only 0.80 cal mol-! K-1. The range 
of existence of the plastic crystalline phase can vary 
considerably. Thus, in the examples given in Table 1, 
it is only about 10º for WF, , but over 200º for camphor. 


A feature of plastic crystals is that they give 
rather poor, indistinct X—ray diffraction patterns, but 
these when analyzed usually reveal that the plastic 
phase is one of high symmetry. Thus, the plastic crys- 
talline phase of cyclobutane only gave two usable lines 
in its X—ray powder photograph, just sufficient to show 
that its lattice is body-centred cubic.39 The low entro- 
pies of fusion and the rather blurred X-—ray diffraction 
patterns are evidence that the crystals are highly dis- 
ordered, and this is amply confirmed by NMR studies. 
An example of the line narrowing associated with a 
plastic crystal has already been given in Fig. 7 for 
pentaerythritol. The line-width just below the melting- 
-point of a plastic crystal is often almost as small as that 
for a liquid. In general, it may be said that in a plastic 
crystal the molecules have considerable orientational 
disorder, and also the ability to move easily from one 
lattice site to another. This ease of movement probably 
cepends on the fact that the crystal has developed an 
appreciable fraction of unoccupied sites (Schottky holes), 
so that a molecule can relatively easily leave its site 
and move into a hole, leaving a hole behind it which 
another molecule can use, and so on. In principle, the 
fraction of such holes can be estimated by comparing 
the actual measured censity of the crystal with that 
calculated from the lattice spacing given by an X-ray 
study, since the presence of holes will mean that the 
actual density will be less than the calculated density. 
It has been found in this way that the percentage of 
unoccupied sites in a plastic crystal does not usually 
exceed one per cent,40 but an exception is provided by 
the plastic crystalline phase of the hydrocarbon cam- 
phene, in which it has been estimated that about seven 
per cent of the lattice sites are unoccupied at the mel- 
ting-point.41 

The adjective 'plastic' reflects the fact that many 
such crystals can be extruded through holes by rela- 
tively low pressures, (e. g. pressures of a few hundred 
atmospheres).38 The plastic crystalline phase of per- 
fluorocyclohexane, C, F,, , is said to flow slowly un- 
der gravity, like pitch. The plasticity of these peculiar 
crystals may well depend, like the ease of diffusion 
of the molecules, on the Schottky holes in the lattice. 

It may be noted that the well-known use of cam- 
phor as a cryoscopic solvent for molecular weight deter- 
minations depends on its existence as a plastic crystal. 
Since the entropy of fusion and the heat of fusion 4 H 
are relatively low, the cryoscopic constant (which is 
inversely proportional to A H f ) is large, and a relati- 
vely large depression of the freezing-point of camphor 
is produced by a relatively low concentration of a solute. 
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G. DISORDER OF POSITION 


In Section A(1), ordinary ice (that is, the hexa- 
gonal form stable at ordinary pressures) was mentioned 
as an example of positional disorder. In addition to 
this stable form, there is an unstable cubic form (also 
disordered with regard to the protons), and no less 
than eight high-pressure forms. Of these, it appears that 
three have the protons in an ordered arrangement, while 
the rest are disordered. The three ordered forms are Il, 
and forms VIll and IX which have only comparatively 
recently been discovered. One simple thermodynamic 
piece of evidence?2 that Il is an ordered structure may 
be mentioned here. From the pressure — temperature 
phase diagram for the ice system (our knowledge of 
which is largely due to the work of Bridgman, supple- 
mented later by the work of Whalley and his associa- 
tes),43 estimates can be made of dp/dT, the slope 
of the line for the co-existence of two phases, such as 
ll and V. The volume change A V for the conversion 
Hi —» V can be measured directly, so from dp/dT, A V 
and the Clapeyron equation 


dp/dT =AS/ AV 


a value is obtained for 4 S for the process Ill->V. 
lt is found in this way that for all possible conversions 
of one phase into another for the forms of ice from | 
to VlIl, those conversions which do not involve Il have 
only small entropy changes, averaging numerically about 
0.07 cal mol-1 K-1, But when Il is converted into any 
one of the three phases which are its neighbours on the 
phase diagram, namely |, Ill and V, (IV is a metastable 
phase), x S is much larger, being (in cal mol! K-1) 
0.8 for ll -- 1, 1.0 for ll -—» Ill, and 1.0 for | — MV. 
Since the entropy which ordinary disordered ice-l 
retains at OK is 0.80 cal mol-1 K-1, the inference is clearly 
that Il is an ordered state, while Ill and V, like |, have a 
disordered arrangement of protons. This conclusion 
about Il is confirmed by its infra-red spectrum (men- 
tioned in Section C(5)), and by other properties such 
as its dielectric constant.44 

Since ordinary ice is disordered at all temperatures, 
it is natural that studies should have been made of 
crystals containing water of crystallization to see if they 
are similar. Experimentally, this is a matter of using the 
procedure outlined in Section B(2), that is of comparing 
the apparent calorimetric entropy of the hydrate, 
S cal with San, the value obtained by an equilibrium 
study. Thus, for Na SO,, 10 H,0,4º Seg can be obtained 
from the equilibrium: 


Na SO,, 10H.0 (c) = Na, SO, (c) + 10 H,0 (9). 
If the water vapour dissociation pressure at, say 25'C 
is PH.O - then AGº = — RTinpg.o - The standard 
enthalpy increase for the reaction, /) Hº, can be obtai- 
ned from the heat of vaporization of water and the heats 
of solution of the anhydrous Na, SO, and the hydrate, 


(strictly, to give very dilute solutions). A Sº for the 
reaction then follows from the relation 


A =" AM =" AS. 
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Using appropriate heat capacity results, the S Nau So, 
the molar entropy of the anhydrous salt is determi- 
ned from the integral |/ Cp /T.dT, (since there 
is no reason to suppose that the anhydrous salt will be 
disordered in any way at OK), while SH,0 , the molar 
entropy of water vapour at one atmosphere and 25ºC, 
is accurately known. San: the true entropy of the de- 
cahydrate, then follows from the equation: 


A Sº — SNa,SO, + 10 SA.o ao “eq 


Seq was found in this way to be 141.46 cal mol! K-1, 
whereas Scal' the apparent calorimetric entropy, is only 
139.95 cal moll K7. So Ss is 1.51 cal moll K'71, 
and hence the hydrate remains disordered at OK. It 
will be noted that S,es is only about one per cent of 
Seq and Seal so that very accurate work is essential. 
A crystal structure study of the decahydrate48 suggests 
that some of the protons of the water molecules in the 
lattice have definite positions, but that for every mole 
of salt there is a group of four water molecules arranged 
in a ring, and the fact that S,es is approximately R In2 
suggests that there are two ways in which the protons 
in the ring can be arranged, as shown in Fig. 12. 
Two other disordered hydrates are Na.CO,,10H 047 
(Sjes — 1.50 cal mol71 K-1) and Na. HPO,,12H.,0,48 
(S,as = 3.51 cal mol-1 K-1). But most hydrates which 
have been studied in this way (mostly by Giauque 
and his collaborators) prove to be ordered. They include 
the following: five hydrates of H SO,; three hydrates 
of NH, ; two hydrates of HNO, ; the hydrates of 
NaOH; ZnSO,,7H,0; Na,CO,,H,O; Na,HPQO,,2H,0 
and 7H, O. 

An interesting and much studied hydrogen-bonded 
crystal is that of KH,PO,. This exists in an ordered, 
orthorhombic form Il and a high-temperature disordered 


form |, the two being separated by a >.-type transition 


approximately tetrahedrally disposed. Similarly, the anion 


HO, /O0 
in KH, PO, , namely P 
HO/ 10 


, has two protons, and 


two oxygen atoms to attract protons, again approximately 
tetrahedrally disposed. Whereas, however, in ice, the 
disordered hydrogen-bonded structure persists to OK, 
in KH, PO, it gives place to an ordered arrangement 
below the transition. The entropy lost when | transforms 
into Il on cooling is 0.74 cal mol-1 K—1. 49 which within 
experimental error is the same as Sres for ice. It then 
seems that as a consequence of the ordering of the 
hydrogen bonds, slight changes in position of the other 
atoms (e. g. the potassium ions) take place, as a result 
of which a unit cell develops a dipole moment. Since 
the structure is ordered, the dipole moment of the cells 
of a group or domain all have the same direction, so 
that the domain of cells is polar. This confers on the 
ordered crystal an exceptionally high dielectric constant 
( - 104); in other words, it belongs to the interesting 
and important class of solids known as ferroelectrics.90 

An example of positional disorder which can be 
described as 'random walk disorder” has recently been 
investigated by Mr. R. D. Worswick in the author's 
laboratory. This occurs in the substance mercury (Il) 
amidobromide, HgNH, Br, which can be obtained in 
two forms, an orthorhombic form and a cubic form.51 
The cubic form is unstable at room temperature, but the 
rate of change into the orthorhombic form appears to 
be so slow that no transformation can be observed over 
a period of months. The cations in both forms consist of 


NH, NH, 
chains Hg Hg Hg in which the 


NH, NH, 


(O O 
De gq 
2 0a, 

O (O 


Fig. 12. Showing the two possibilities for hydrogen-bonding in a ring of four water molecules in the 
lattice of Na,SO,, 10H O. The larger circles are oxygen atoms, the smaller circles are hydrogen atoms. 


with its heat capacity maximum at = 123K. There is a 
similarity here, perhaps at first sight unexpected, with 
ice-l. The structure of ice is determined by the fact 
that each water molecule has two protons for forming 
hydrogen bonds and two lone pairs of electrons for 
attracting the protons of other molecules, the two pro- 
tons and the two lone pairs of a given molecule being 
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unit N—Hg—N is linear. In the orthorhombic crystal, the 
Hg and N atoms of any cation chain lie in a plane, for- 
ming a regular zig-zag, and adjacent chains are parallel 
so that all Hg atoms have definite positions. In the 
cubic form, the N atoms lie at the cube corners, and 
every Hg atom must be situated half-way between two N 
atoms, that is at the mid-points of the cube edges. 
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Only one-third of these possible sites for the Hg atoms 
is occupied, however, so that the Hg atoms are posi- 
tionally disordered. This is equivalent to saying that the 
path taken by a cation chain in the cubic form is a path 
followed on a walk through the cubic lattice in which 
every step (i. e. an N—-Hg-—N unit) is an edge of a 
cube, and in which each step must be at 90º to the 
preceding step. The determination of the number of 
such walks which are possible in a given crystal is 
another example of a configurational problem which 
cannot be exactly solved, but an approximate estimate 
has been made by Dr. D. F. Mayers by computing the 
number of random walks possible for blocks of unit 
celis of increasing size, beginning with 3x3x3. The 
computation has been carried as far as blocks of 
13x13x13 unit cells, (the computing time required rising 
rapidly, of course, with increasing size of the block). 
Extrapolation of the results so obtained to (in effect) 
the infinitely large crystal, gives an estimate of about 
1.7 cal mol K-1 for the entropy associated with 'ran- 
dom walk disorder" in a mole of Hg NH, Br. The ex- 
perimental examination of whether or not the crystal 
retains entropy of about this amount at OK has been 
made by comparing its apparent calorimetric entropy 
Scal at 25 C with Se Seq was determined by measu- 
ring the molar entropy at 25ºC of the orthorhombic 
form, Sor' and also the difference Sor Seq: This 
difference was cetermined with the help of the cell 
Hg | HoNH, Br (c) | KBr (aq) | HgNH, Br (c) | Hg, 
(cubic) (orthorh.) 
in which the cell “reaction” is 


HgNH, Br (cubic) = HgNH, Br (orthorhombic). 


This cell gave a fairly reproducible E. M, F. of 0.0050 
volts, which corresponds to 4 Gº = — 230 cal mol-?. 
/ Hº for the cell reaction was found to be —15.5 cal 
mol-1 from the difference in the heats of solution of the 
two crystalline forms in hydrochloric acid. These results 
give sa Sº = Sor me Seq = 0.72 cal mol-1 K-1. Scal for 
the cubic and orthorhombic forms were determined by 
appropriate low-temperature calorimetry to be 31.31 and 
31.83 cal mol-! K-1 respectively. Sos for the cubic 
form then follows from the equation 


s + 31.31 + 0.72 = 31.83, 
res 


giving Sjes = — 0.20 cal mol- K-1, with an estimated 
uncertainty of + 0.3 cal mol-! K-1. Therefore S,as is 
zero within experimental error. This is a somewhat sur- 
prising result, since it means that the cubic form, while 
undoubtedly disordered at room temperature, somehow 
contrives to become ordered (as regards the position 
of the Hg atoms) on cooling, although there was no 
sign of any transition in the heat capacity measure- 
ments, which extended from - 10K to 300K. 


H. MAGNETIC DISORDER 


This is a large subject, and we shall limit ourselves 
to a brief illustration of a point not so far made, 
namely that the order-promoting influence in a crystal 
may be much stronger in one direction than in an- 
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other. Consider, for example, a copper (Il) compound. 
The ion Cu2+ with a 3dº structure has one unpaired 
electron, so that the degeneracy which can result from 


1 
the spin of S = EA is 2S + 1, = 2, and the total ma- 


gnetic entropy gained by the spin system becoming 
disordered is therefore R In2. The ordering at sufficiently 
low temperatures can be antiferromagnetic or ferro- 
magnetic. In Cu Cl,, the transition from the antiferro- 
magnetically ordered state to the disordered state gives 
a pointed maximum in the heat capacity at about 24K.52 
At this temperature, however, only about 27 per cent 
of the magnetic entropy of R In2 has been acquired. 
Analysis of the heat capacity shows that the rest of 
the magnetic entropy is developed between 24K and 
about 100K. This is readily interpreted in terms of the 
crystal structure of CuCl, , which consists essentially 
of parallel chains: 


Cl Cl Cl CI 
Cu Cu Cu | | Cu 
Cl Cl Cl Cl 


The magnetic interaction between the copper atoms 
along the chain is much stronger than that between 
copper atoms of different chains. At sufficiently low 
temperatures the system is magnetically ordered in 
three dimensions, but the magnetic order between the 
chains is relatively easily broken down, while leaving 
considerable order along the chains which finally dis- 
appears only at much higher temperatures. This is no 
doubt a case where the interaction between the copper 
atoms, particularly along the chains, depends on super- 
exchange (Section D). 

We have recently investigated magnetic transitions 


in double copper chlorides containing the anion 
Cu, Cl, , Of structure 
Cl Cl Cl 
Cu Cu 
| ai 
Cl Cl Cl 


An example of such a salt is Li Cu Cl,, 2 H,0, which 
undergoes a gradual antiferromagnetic -» paramagnetic 
transition with a maximum in the heat capacity at about 
4K.53 The entropy gain associated whith this transi- 
tion was found to be 1.07 + 0.02 cal mol-1 K-1, which 


is less than R In2, but significantly, almost equal to 


: R In3. The interpretation of this result is that the two 


spins of the copper atoms within one and the same 
anion are strongly coupled in a parallel way, giving a 
resultant spin S = 1, and a degeneracy of 3. There 
are N/2 anions per mol of the double chloride, so if 
it is supposed that above 4K the magnetic order between 
the anions has disappeared (while the coupling between 
the two copper atoms in any one anion still persists), 
the expected magnetic entropy gain shoul be k In 3N/2, 


É 
ad R In3, in agreement with the experimental result, 
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Thus the conclusion reached from this simple thermo- 
dynamic observation is that in the Cu.Cl ion, the 
superexchange is of a kind which produces a relatively 
strong interaction between the copper atoms, resulting 
in parallel spins of the unpaired electrons. 


Permission to reproduce the following figures is 
gratefully acknowledged: fig 3 (from ref. 9) by permis- 
sion of the Chemical Society of Japan; figs. 4, 7 and 
9 (from refs. 12, 17 and 21 respectively) by permission 
of the authors; fig. 5 (from ref. 14) by permission of 
the Editor of 'Nature' and the authors; fig. 6 (from 
ref. 15) by permission of the author, the Institute of 
Physics and the Physical Society; fig. 8 (from ref. 20) 
by permission of the North-Holland Publishing Com- 
pany; fig. 11 (from ref. 34) by permission of the Ame- 
rican Chemical Society. 
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LIFTING LINE THEORY OF AXIAL FLOW TURBOMACHINES — À REVIEW (*) 


SYNOPSIS 


The paper presents a review of the literature on 
the use of lifting line theory in the analysis of the flow 
through axial turbomachines. The case of screw propellers 
and open turbomachines, the effects of cylindrical walls, 
compressibility and rotational flow are successively 
considered. 


1. — INTRODUCTION 


In the design of axial turbomachines, the analysis 
of the flow through a blade row is frequently carried 
out as the combination of two different approaches: 

a) Blade element flow, which enables the designer 
to use the performance data obtained mainly from two- 
dimensional flow investigations with isolated aerofoils 
or with linear cascades. It is assumed that the flow 
about any blade section remains the same regardless of 
whether the section is in a two-dimensional plane flow 
field or in an arbitrary annular configuration. 

b) Flow induced by the whole blade row or rows. 
Due to difficulties of the full three-dimensional flow 
field analysis, simplifying assumptions are always made, 
leading to different theories. Axial symmetry is one of 
the most important assumptions, by which the partial 
derivatives of the flow variables with respect to the 
angular coordinate are taken equal to zero. Physically, 
this implies that the blade row is constituted by an 
infinite number of infinitely thin blades and the shed 
vorticity (which occurs whenever the circulation around 
the blade section varies in the spanwise direction) is 
circumferentially uniform. 

There are effects which are not accounted for in the 
combination of axially symmetric theory and two-dimensi- 
onal blade element theory. One of these is the fact that 
the trailing vorticity from the actual blade row is concen- 
trated in the blades' wakes rather than circumferentially 
uniform. It is expected that this modification to the 
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RESUMO 


Apresenta-se uma revisão bibliográfica sobre a 
aplicação da teoria da linha sustentadora ao estudo do 
escoamento em turbomáquinas axiais. São sucessiva- 
mente considerados o caso dos hélices propulsores e 
outros tipos de turbomáquinas abertas e os efeitos devi- 
dos à presença de paredes cilíndricas, compressibilidade 
e rotacionalidade do escoamento. 


induced velocity field should be particularly important 
for relatively small numbers of blades and for conside- 
rable variation of the circulation along the span. 


2. — LIFTING LINE THEORY 


The lifting line theory of axial flow turbomachines 
follows from Prandtl's theory of wings of finite span 
[1]**. The blades are replaced by bound vortex lines 
or lifting lines, which are stationary in the case of linear 
cascade or stator blades, or are supposed to rotate 
about a fixed axis in the case of rotor blades, and the 
velocity circulation at any point of the lifting line is set 
equal to the circulation about the corresponding blade 
section. 

This simplification is expected to yield good appro- 
ximations for the resultant forces exerted by the flow 
upon the blades and to give an approximate picture of the 
overall flow induced by the blade row,although the 
flow in the vicinity of a lifting line is definitely different 
from the flow near a blade surface. In particular, the 
method can provide a good tool for the study of the 
interference effects due to the blades and their wakes, 
but accurate results should not be expected if the blade 
chord is not small compared with the distance between 
neighbouring blades. 

When the velocity circulation — or the blade load — 
is not constant along the span, there must be, according 
to Kelvin's theorem, a wake constituted by free vortex 


(* *) For more recent accounts of Prandtl's theory, see, for example, [2 to 4]. 
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filaments shed by the blade. In the simplified lifting 
line model, the wakes take the shape of vortex sheets. 
For isentropic flow, the second vorticity theorem of 
Helmholz (see e. g. [5]) states that the trailing vortex 
filaments move with the fluid and, if the flow is steady, 
coincide with the streamlines. 

It follows that the motion of the shed vortex fila- 
ments and their instantaneous shapes depend upon the 
velocity field they induce themselves. The exact treat- 
ment of this case leads to a non-linear mathematical 
problem, the solution of which involves great analytical 
and numerical cifficulties, even for the trivial case of a 
statiorary straight lifting line with a trailing vortex 
sheet in an undisturbed uniform flow. 

in the classical approach to the problem, it is 
assumed that the trailing vortex filaments have simple 
geometrical configurations (straight or helical lines), 
usually coinciding with the primary flow streamlines 
(or the relative streamlines, in the case of a rotor). This 
assumption implies that the disturbance — or velocity 
induced by the bound and trailing vortex systems — is 
small compared with the undisturbed flow velocity, so 
that the actual streamlines do not differ considerably 
from the streamlines of the primary flow. We obtain, in 
this way, the so called first order flow based on small 
disturbance linear theory. Prandtl's lifting line theory of 
wings with finite span is the best known example of 
this approach. 


3. — SCREW PROPELLERS AND OPEN TURBOMA- 
CHINES 


It is not surprising that the method was first applied 
to screw propellers, where the circulation, being zero 
at the blade tip, has a great spanwise gradient and the 
shed vorticity, concentrated in helical wakes, induces 
a velocity field about the blades which cannot be negle- 
cted without introducing considerable errors. Moreover, 
as the number of blades is usually small, actuator disc 
or any other theory based on axial symmetry can only 
provide a rough approximation to the actual flow. 

Three types of approximation are usually considered 
in propeller vortex theory. In the case of a lightly loaded 
propeller, the influence of the velocity induced by the 
bound and trailing vortex systems on the shape of the 
vortex sheets is neglected. The free vortices are then 
helices and the vortex sheets constitute true helical 
surfaces, the pitch of which is determined solely by 
the ratio of the advance speed to the angular speed. In 
the case of a moderately loaded propeller, the free vortex 
lines are still suposed to lie on cylindrical surfaces of 
revolution but the influence of the tangential and axial 
components of the induced velocity upon the pitch is 
now taken into account and the vortex sheets are not, in 
general, true helical surfaces. Finally, for a heavily 
loaded propeller, the effect of the centrifugal forces 
and the contraction of the slipstream are also taken 
into account. 

Betz [6], in 1919, applied, to screw propellers, 
Prandtl's lifting line theory of wings of finite span and 
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established the condition for minimum energy loss at a 
given thrust for a moderately loaded propeller: the flow 
far behind the propeller must be such that the helical 
surface swept by each blade moves downstream like 
a rigid body, i. e. the pitch of the trailing vortices does 
not vary with the distances to the axis of the propeller. 
By applying certain symmetric properties of the vortex 
system, Betz showed that the downwash velocity (obtai- 
ned by combining the tangential and axial componets of 
the velocity) induced at the spoke lifting lines is one 
half of the velocity induced at corresponding points 
of the vortex sheets far downstream. This important 
relation was later widely used by other research workers 
and designers. In an appendix to the same paper, Prandtl, 
replacirg the helical vortex sheets by parallel plane 
sheets, derived an approximate solution for the optimum 
radial distribution of the circulation. 

In 1929, Goldstein [7], using Betz's conclusions, 
solved Laplace's differential equation for the velocity 
potential far downstream in a non-orthogonal helical 
system of coordinates and obtained the radial distribution 
of the circulation which gives minimum energy loss. In 
the same paper, numerical results were presented for 
a two-bladed propeller. More complete evaluations for 
several numbers of blades were later given by Kramer 
[8] and, more recently, by Tatchmindji and Milam [9]. 

Kawada [10, 11], using Goldstein's helical system 
of coordinates, solved the potential problem far behind 
a screw propeller with N blades, i. e. for a system 
constituted by N lifting lines of constant circulation, 
N free helical vortex lines shed by the blades' tips and 
an axial straight vortex line of the combined strenght 
of the lifting lines, extending downstream from the centre 
of the propeller. The solution involved series in terms 
of the modified Bessel functions | and K, which, for 
numerical evaluations, were replaced by Nicholson's 
asymptotic expressions. 

The solution of a similar mathematical problem had 
been published before by Lamb [12] for the magnetic 
field of a helix. 

Moriya [13], by direct integration of Biot-Savart's 
equation and also under the assumption of constant blade 
circulation, obtained the expressions already derived by 
Kawaca for the downwash velocity at points of a lifting 
line. 

Once the velocity field induced by a free helical 
vortex line of radius É is known — this is the problem 
solved by Kawada — the field, due to a lifting line of 
variable circulation and its trailing vortex sheet, can be 
obtained by integration with respect to : between * =r, 
and ro (respectively the radial coordinates of the blade 
root and tip). However, a complication arises when the 
downwash velocity at a point P(r) of the lifting line is 
to be calculated by this method. In fact, as the radius 
of the helical line * approaches the radial coordinate r 
of P — i. e. as the starting point of the helical vortex 
tends to P — the velocity induced at P tends to infinity, 
which means that the integrand has a singularity at? = r 
and the integral exists only as a Cauchy principal value. 
This fact obviously makes difficult the application of qua- 
drature methods. 
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Moriya [14] suggested a method for evaluating the 
integrais which give the tangential and axial components 
of the velocity induced at a lifting line of variable circula- 
tion, by introducing the so called tangential and axial 
induction factores. The tangential induction factor is 
defined by the ratio 


v, (1, E.B) 


wir,5) 


| (r,E.B) = 


where v, (r,E,8) is the tangential component of the 
velocity induced at point P(r) by a free helical vortex 
of unit strength, with radius É and angle É. In the 
denominator, w (r,;) is the velocity (with tangential 
direction) induced at P (r) if the helical line is replaced 
by a semi-infinite straight vortex of unit strength, para- 
lel to the axis of the propeller and extending to infinity 
from the same starting point. Biot-Savart's law gives 


| 1 
wir, ES RS 
(r, E) pp 
It is easy to see that i, (r, É, 6) remains finite, even 
foré = r. 
The tangential component of the velocity due to 
the whole vortex sheet is then given by 


1 fe dV'(E) i, E E 8) 


Vo (r) = 
4r q F; dE , — 

where |" (?)is the circulation about the lifting line and 8 
is assumed to be a known function of É. The right hand 
side of the last equation is formally identical to the ex- 
pression of the downwash velocity of a wing with a span 
equal to rg — r;, the equivalent spanwise distribution 
of the circulation ['*(:) being such tha d['*/d£ 
=, ql dt 

A similar analysis applies to the axial component 
of the induced velocity with the aid of the axial induction 
factor ia. 

Provided that the induction factors are known, for 
each value of r, as functions of :, the integration can 
be carried out with the help of methods derived in the 
theory of wings of finite span. Lerbs [15, 16], extending 
Glauert's aerofoil method [3], determined the principal 
value of the integrais by expanding the circulation |" 
and the induction factors in Fourier series in terms of a 
variable related to E. 

Values of the induction factors were published in 
the form of diagrams by Moriya, who evaluated, by a 
graphical method, the integrals arising from the applica- 
tion of Biot-Savart's law, for 2-bladed propellers [14] 
and later also for 3- and 4-bladed propellers [171]. 
More detailed and accurate diagrams for 3-bladed pro- 
pellers were later published by Strscheletzky [18, 19], 
who obtained the results by direct numerical integration 
of the Biot-Savart's integrals. Values of the induction 
factors, for propellers with 3, 4 and 5 blades, were 
obtained and published in the form of diagrams by Lerbs 
[15, 16]. who used Kawada's analytical solution, with 
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Nicholson's asymptotic expressions replacing the modi- 
fied Bessel functions | and K. More accurate and exten- 
sive numerical results, presented in tables, were later 
obtained by Morgan [20, 21] with the aid of more 
approximate asymptotic expressions of the modified 
Bessel functions. 

A different way of dealing with Biot-Savart's inte- 
grals was used by Guilloton [22, 23], who wrote the 
velocity induced at a point r of a lifting line as 


1 fr dD6 d 1 Pro dT() 
Wipe] —— + — ) Hai 
4xJ" di E-r 4x fo di 


where H is a function of r, É and the helical pitch, 
which represents the effect of the helical shape of the 
trailing vortices (in comparison with the straight vortices 
assumed in wing theory) and remains finite when £ —»r. 
Only the first integrand presents a singularity at £ = r 
and its integration can be overcome by Glauert's method 
[3] for wings of finite span. In fact, there still remains, 
in H, a singularity of the type of In (£ — r), which was 
found by Guilloton to be practically negligible. Numerical 
values of the second integral were given by Guilloton 
in terms of the first two coefficients of the Fourier de- 
velopment of I'(:), for several combinations of the 
geometrical parameters of the propeller. 

In all the works refered to, where the concept of 
induction factor is used, the slipstream contraction and 
the effects of the centrifugal force are neglected, but the 
variation of the pitch of the helices with the radius can 
be taken into account. This means that those methods 
apply to moderately or to lightly loaded propellers. 

A comparison of the results obtained by using 
some of these lifting line methods in propeller design 
was presented by Johnsson [24]. 

A different approach was used by Reissner [25, 26], 
who solved the potential problem in the far wake of a 
screw propeller for an arbitrary spanwise distribution of 
the circulation, assuming that the trailling vortex sheets 
constitute true helical surfaces. To obtain the appropriate 
discontinuities of the potential and of the radial com- 
ponent of the velocity across the vortex sheets, he used 
an artifice which consists in introducing a source dis- 
tribution of helical symmetry, with a linear saw-tooth 
circumferential variation, which produces the desired 
discontinuities. A second source distribution is therea- 
fter superposed, in such a way that its circumferential 
variation represents a Fourier series of the saw-tooth 
curve, Poisson's differential equation is solved for both 
distributions, in a system of helical coordinates identical 
to that used by Goldstein and Kawada, and, subtracting 
the second solution from the first one, the flow field 
induced by a set of equidistantly spaced helical vortex 
sheets is obtained. Having found the expressions for the 
downwash velocity induced at the lifting lines (which 
is equal to one half of the induced velocity at correspon- 
ding points of the vortex sheets far downstream, accor- 
ding to Betz's law for semi-infinite helical vortex sheets), 
Reissner combined these results with blade element two- 
dimensional theory to establish a relationship between 
the blade geometry and the induced downwash velocity. 


193 


Further developments of the vortex theory of the 
screw propeller in incompressible flow, including lifting 
surface theory, can be found in [27 to 30]. 


4. — THE EFFECT OF CYLINDRICAL WALLS 


To evaluate the velocities induced at the blades 
of a propeller operating in free space, Biot-Savart's 
equation can be applied without great difficulties and 
with good results. However, when adjacent boundaries, 
such as a hub, a duct or both, are present, the condition 
of zero normal velocity at the walls must be satisfied 
and the use of Laplace's or Poisson's differential equation 
for the velocity potential is a more suitable approach, par- 
ticularly if the walls are cylindrical surfaces. Obviously, 
in this case, the analytical solutions become more com- 
plex. 

Takeyama [31] extended Kawada's analytical 
solution to the case of N symmetrically spaced helical 
vortices bounded by a circular cylindrical casing exten- 
ding to infinity in both axial directions. Additional terms 
to Kawada's expressions were derived by Lerbs [15] to 
take into account the effects of a cylindrical hub of 
infinite length. 

Hickling [33] presented the expression of the 
velocity potential for the case of N symmetrically spaced 
helical vortices enclosed in an annular cylindrical duct 
of infinite length and also gave the simplified expressions 
for the cases when only one of the walls is present. 
The complete solution is given by an infinite series 
with terms involving cross products of the modified 
Bessel functions | and K and their first derivatives. With 
the help of numerical examples, Hickling concluded that 
the flow field in the midde of the annular duct is only 
slightly affected by the presence of the walls and that, 
near the hub or the casing, the solution approaches that 
of the case when only one boundary is present. With 
an increase in the number of blades, the effects of the 
hub and duct walls become confined to smaller regions 
and, for an infinite number of blades, the boundaries 
do not affect the flow. 

Meyer [34], in 1946, studied the flow in a statio- 
nary annular cascade of radial source lines and lifting 
lines of constant circulation, by integrating Laplace's 
potential equation in cylindrical polar coordinates by the 
method of separation of variables. The solutions are 
expressed by double series involving Bessel functions J 
and Y of the radial coordinate and circular functions of 
the angular coordinate. The effect due to the finite num- 
ber of blades in shown to decay exponentially away from 
the cascade plane. Meyer also extended, to a stationary 
annular cascade of blades with constant circulation, 
Ackeret's [35] singularities method for two-dimensional 
linear cascades. 

The deviation of the two-dimensional approxima- 
tion (i. e. linear cascade) from the three-dimensional 
annular cascade of lifting lines of constant strength was 
studied theoretically by Ginsburg [36], who considered 
particularly the flow field in the vicinity of the cascade. 
Using an electrolytic tank, he measured the potential 
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field for a hub/casing diameter ratio of 0,308 and 
compared the results with the theoretical values. 

Rossow [37] extended Ginsburg's theoretical in- 
vestigation to the region away from the cascade for the 
cases of radial source lines and radial lifting lines of 
constant strength. In his analysis, the three-dimensional 
solution is expanded in a Taylor series about the limi- 
ting case of hub/casing diameter ratio equal to ore. 

The effects of the trailing vorticity in an annular 
cascade were studied by Nath, who carried out the 
numerical integration of Biot-Savart's equation for three 
cdifíerent cases: straight axial vortex filaments concen- 
trated along symmetrically spaced plane sheets [38]; 
helical vortex filaments of constant pitch axisymmetrically 
cistributed between coaxial cylinders and shed by an 
actuator disc [39]; and finally the case of helical vortex 
sheets equidistantly spaced [40]. The circulation about 
the lifting lines was taken equal to a third degree poly- 
nomial of the radial coordinate, with the first derivative 
vanishing at the inner and outer walls (in order not to 
have infinite downwash velocities at the ends of the 
lifting lines). However, in his three papers, Nath, although 
considering the problem of ducted fans, neglected the 
effect of the inner and outer boundaries and, as a conse- 
quence, strictly his results (with the exception of the 
actuator disc) apply only to a propeller or a fan in free 
space. 


5. — COMPRESSIBLE FLOW 


In all the works quoted so far, only the case of 
incompressible flow was considered. The vortex theory 
of the screw propeller in compressible flow was sketched 
by Busemann [41]. Davidson [42], using Reissner's 
analysis [25, 26] for the potential in the far wake, 
presented a solution for the whole flow field for a 
propeller of arbitrary circulation in compressible flow. In 
Davidson's analysis, the propeller is enclosed in a 
coaxial cylindrical wind tunnel of circular cross section, 
supposed to extend to infinity in both axial directions. 
It was later pointed out by Okurounmu [43] that his 
choice of the boundary conditions at the cascade plane 
is inappropriate and the solution does not satisfy the 
continuity requirements for compressible flow. 

Following the analysis of McCune [44, 45] for 
the three-dimensional compressible flow in an annular 
cascade of source lines, Okurounmu and that author 
[43, 46] studied the case of a rotating annular cascade 
of lifting lines of prescribed variable circulation in 
compressible flow. Using a non-orthogonal helical system 
of coordinates, as Goldstein, Kawada and Reissner had 
done for the screw propeller, they showed that, far 
downstream from the rotor, the velocity field can be 
expressed in terms of only two coordinates — the 
radius and the helical angular coordinate — not explici- 
tiy dependent on the axial coordinate, and the solution 
has the same form as for the case of incompressible 
flow. The expressions for the potential far downstream 
were obtained by extending Reissner's results to take 
into account the presence of cylindrical walls, and the 
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expressions for the complete flow field were then obtai- 
ned by superposing, on the far wake solution, a true 
three-dimensional solution that decays exponentially to 
zero with the distance to the cascade plane and which 
includes the compressibility effects. 


6. — ROTATIONAL FLOW 


in all the works refered to before, either those 
using Biot-Savart's law or based on the velocity potential, 
it is explicitely or implicitely assumed that the distur- 
bance flow field is irrotational everywhere except at the 
singularities (the bound and trailing vortex lines). As we 
show below, this assumption implies that the primary or 
undisturbed flow — i. e. the flow in the absence of 
blades — is irrotational. However, in many problems of 
fluid mechanics in turbomachines, the actual flow field 
differs considerably from the conditions of irrotationa- 
lity, due to the presence of wall boundary layers or to 
the passage through stationary or rotating rows of blades. 

If the primary flow is not irrotational everywhere, 
the effect of the disturbance due to the blades is to alter 
the pattern of the streamlines and, consequently, also 
the pattern of the primary flow vortex lines convected by 
the flow. 

If the velocity disturbance is of order (=) (the 
primary flow velocity Vo being of order of unity and s 
being a small quantity), the disturbance to the vorticity 


vector field is of order of «(Jo (where Do = curl Vo 
is the vorticity vector in the primary flow). It follows 
that, unless it is Oo = O everywhere, the disturbance 
flow outside the bound and trailing vortex systems is 
rotational, with vorticity of order of «Oo. 

lf the primary flow vorticity Oo is everywhere 
small, say O(:), we have s0o = O(s?) and, negle- 
cting quantities of O( :*) or smaller, we can take the 
disturbance as irrotational. This is the first approxima- 
tion in Hawthorne's theory of shear flow [47, 48] — 
small shear, small disturbance. 

If, however,Qo is O(1), the vorticity of the dis- 
turbance flow is O( :), i.e., is of the order of the velo- 
city disturbance itself and cannot be neglected(*). The 
secondary flow field must, therefore, be considered as 
rotational. This case corresponds to Hawthorne's second 
approximation — large shear, small disturbance. 

Some solutions concerning small perturbation due 
to lifting systems in parallel shear flow can be found in 
the literature. One of the earliest contributions is the 
work by von Kármán and Tsien [49], who extended 
Prandtl's lifting line theory to the more general case of 
a rotational primary flow. Honda [50], using a lifting 
surface theory, obtained a solution for the case of 
incompressible shear flow about a thin wing spanning 
two parallel walls and later [51] extended the solution 


(*)The vorticity being given by a first order differential 
operator of the velocity, it follows that a vorticity field of Ota) 
will induce a velocity field of the same order of magnitude if 
the representativo spatial dimensions of the field (e.g., the 
cascade pitch) are O(1). 
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to the case of a linear cascade of thin blades in incom- 
pressible shear flow. Numerical results for this last pro- 
blem were obtained by Nally and Hawthorne [52] with 
the aid of a digital computer. Namba and Asanuma 
[53, 54] studied the case of a rectilinear cascade of 
lifting lines in subsonic shear flow and later Namba 
[55] extended Honda's incompressible flow analysis to 
the case of a linear cascade of thin blades in subsonic 
shear flow. A more general and detailed analysis of 
shear flow problems can be found in a work by Hawt- 
horne [471]. 
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SISTEMA PARA MEDIÇÃO DE TEMPOS CURTOS DE 
PROPAGAÇÃO DE ONDAS SÍSMICAS 


RESUMO 


Trata-se dum sistema para medição do tempo de 
percurso duma frente de onda entre dois pontos afas- 
tados de 1 a 5 metros, visando a caracterização dum 
maciço rochoso do ponto de vista de resistência mecá- 
nica. Referem-se em especial as unidades electrónicas 
desenvolvidas no LNEC e comparam-se as suas caracte- 
rísticas com as do equipamento corrente de prospecção 
sísmica. 


1 — INTRODUÇÃO 


A caracterização de formações geológicas do ponto 
de vista de resistência mecânica pode efectuar-se a par 
tir das respectivas velocidades de propagação obtidas 
em ensaios sísmicos. A exploração desta potencialidade 
tem especial interesse quando efectuada em conjunto 
com ensaios estáticos, tendo em vista o estudo com- 
parativo dos resultados. 

Para que os dois métodos sejam correctamente 
correlacionáveis torna-se porém necessário adoptar um 
dispositivo sísmico de dimensões suficientemente pe- 
quenas para que seja respeitado o carácter local do 
ensaio. Desta circunstância resulta um método de me- 


EXPLOSOR 


J. GUIMARÃES CORREIA 
Engenheiro Electrotécnico, Especialista 
do Serviço Técnico Geral do LNEC 


SYNOPSIS 


This work concerns a system for measuring the 
travel time of a wave front between two points distant 
from 1 to 5 meters with a view to characterizing a rock 
mass from the viewpoint of mechanical strength. The 
LNEC's electronic units are described in more detail and 
their characteristics are compared with those of the 
current equipment for seismic exploration. 


dida, parcialmente diferenciado em relação ao da pros- 
pecção sísmica clássica, a que se dá a designação de 
microssísmico e para o qual foi desenvolvido o equipa- 
mento que se apresenta nesta memória. 


2 — TÉCNICA DE MEDIDA DO MÉTODO MICROS- 
SÍSMICO 


Na figura 1, apresenta-se o esquema de medida do 
método microssísmico, A actuação do explosor faz deto- 
nar uma cápsula, originando uma onda sísmica que é 
captada sucessivamente por dois geofones. Os respec- 
tivos sinais eléctricos, depois de amplificados, são exi- 
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Fig. 1 — Esquema de medida do método microssísmico 
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bidos nos dois canais do osciloscópio, cuja base de 
tempo é disparada, com ou sem atraso, a partir dum 
impulso fornecido pelo explosor. Há apenas uma passa- 
gem do feixe que é fixada fotográficamente de modo a 
permitir avaliar o tempo de percurso da frente de onda 
entre os dois pontos de captação. 

Explosor e geofones são aparelhos que fazem parte 
do equipamento corrente de prospecção sísmica. 

Em relação a este a principal diferença reside na 
unidade de registo. A necessidade da sua substituição 
pelo sistema osciloscópio câmara fotográfica resulta de 
serem bastante mais curtos os intervalos ce tempo a me- 
dir. Com efeito, a distância entre geofones poderá variar 
entre 1 e 5 metros, a que, para uma velocidade de, por 
exemplo, 2000 metros por segundo, correspondem tempos 
de percurso entre 0,5 e 2,5 milissegundos. Para se medir 
com suficiente resolução dentro desta gama é necessá- 
rio um sistema de registo com uma base de tempo bas- 
tante mais rápida do que podem proporcionar os apare- 
lhos de prospecção sísmica. 

Quanto aos amplificadores e retardador de impul- 
sos, dado que se construíram unidades especialmente 
destinadas ao método microssísmico, serão objecto dum 
tratamento mais pormenorizado. 


3 — AMPLIFICADORES SÍSMICOS 
3.1 — Sobre as características do sinal microssísmico 


No método microssísmico poderiam, em princípio, 
utilizar-se amplificadores sísmicos clássicos. Esta solu- 
ção não é contudo a mais conveniente, pois existem 
certas diferenças quanto às caracteristicas mais apropria- 
das para cada um dos tipos de ensaio. 

No que respeita a ganho, as necessidades do mé- 
todo microssísmico são bastante inferiores aos máximos 
disponíveis nos amplificadores de prospecção. Este facto 
resulta das menores dimensões do dispositivo, pois, em 
regra, a distância do ponto de tiro ao geofone mais 
afastado não ultrapassará a dezena de metros. 

Quanto à resposta em frequência, um amplificador 
de prospecção está limitado inferior e superiormente a 
valores da ordem de 10 ciclos por segundo e 300 
ciclos por segundo, banda que não é a mais conveniente 
para o método microssísmico. 

De facto, o espectro de amplitude dum impulso 
sísmico captado à distância x do ponto de explosão 
pode considerar-se da forma [1]: 


— a (mw) X 
Ax (0) = Ag (mv) 


em que a constante de amortecimento z (| w|)é crescente 
com w e dependente das características do meio, dimi- 
do com o seu grau de consolidação. 

Atendendo a que a prospecção sísmica se aplica 
em geral em meios heterogéneos, incluindo com fre- 
quência um solo, e que o método microssísmico tem em 
vista o ensaio de maciços rochosos, a constante « será 
em regra mais pequena neste caso. Como, por outro 
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lado, as distâncias são também menores, o mesmo amor- 
tecimento ocorre no segundo caso a uma frequência 
mais alta do que no primeiro. 

Conclui-se portanto que o espectro do impulso 
típico de método microssísmico ocupa uma banda de 
largura superior à do sísmico. 


3.2 — Características dos amplificadores 


A anterior conclusão foi confirmada por experiên- 
cias que permitiram especificar, para projecto dos am- 
plificadores, uma largura de banda compatível com uma 
boa relação sinal - ruído. Também por via experimental, 
foi estabelecido um valor conveniente para o ganho. 

O dimensionamento foi ainda condicionado a um 
baixo factor ce ruído, baixa impedância ce saída, con- 
sumo reduzido e uma impedância de entrada relativa- 
mente baixa, adequada às características dos geofones. 
O circuito inclui um atenuador de entrada, dois andares 
amplificadores e um andar de colector comum. 

As características dos amplificadores (unidade 
AS-03, frequência da ensaio 500 c/s) são as seguintes: 


Ganho: máximo 530, em escalões nominais de 500, 
200, 100, 50, 20 e 10. 

Largura de banda: — 3 dba 59 c/s e 6,6 kc/s 
impedância de entrada: 950 º 

Impedância de saída: 140 0 

Amplitude máxima de saída: 5V, pico a pico 
Tensão ce ruído equivalente à entrada (gerador - 
- 220 “9, ganho — 500): 1» W, valor eficaz 

Tensão de alimentação: 9 V 

Consumo: 2 mA 


Os amplificadores foram construídos com a forma 
cilíndrica (figuras 2 e 3) que é apropriada a uma apli- 
cação com muito interesse do método microssísmico: 
o ensaio de maciços em profundidade, ao longo de 


Fig. 2 — Aspecto exterior do amplificador 


furos de sondagem. Esta forma permite a introdução dos 
aparelhos nos furos de modo a encurtar, com vantagem 
do ponto de vista de ruído, os cabos de ligação aos 
geofones. 

Foi também projectada uma blindagem estanque 
(figura 4) que protegerá os amplificadores nos casos, 
bastante frequentes, em que os furos contenham água. 


TÉCNICA 414 


